Going beyond GHZ paradox 
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We present numerical data showing, that three qutrit correlations for a pure state, which is not 
maximally entangled, violate local realism more strongly than three-qubit correlations. The strength 
of violation is measured by the minimal amount of noise that must be admixed to the system so 
that the noisy correlations have a local and realistic model. 

<n : 

The seminal paper of Greenberger, Home, and Zeilingcr |l[ has initiated a completely new phase in the discussions 
regarding the Bell theorem Q . Einstein-Podolsky-Rosen |3| elements of reality were suddenly ridiculed by a straight- 
forward argumentation. The physics community immediately noticed that the increasing complexity of entangled 
systems does not lead to a less pronounced disagreement with the classical views, but just the opposite! Moreover, 
■ the disagreement exponentially grew with the number of qubits involved in the GHZ-type entangled states. Indeed, 
prior to the publication of Ref. [Q , it was commonly perceived that everything regarding the Bell theorem is known. 
C*"} ■ However, the new insight has renewed the interest in the Bell theorem and its implications. 

Another widely shared perspective was that one cannot gain additional useful insight into the Bell Theorem by 
increasing the dimensionality of the entangled systems. Some papers even suggested that in N dimensional systems, 
] increasing the dimension N effectively brings the system closer and closer to the classical realm. However, due to the 
t-H ■ fact that the N > 2 dimensional systems can reveal the Kochen-Specker paradox , this view could be challenged. The 
advent of the quantum information theory created the awareness that such systems require much less entanglement to 



be non-separable than qubits Q . Certain strange features like bound entanglement | or inextensible product bases 



X 



I jjj, suddenly emerged. 

Recently, it was shown that higher dimensional entangled systems indeed may lead to stronger violations of local 
realism, even in straightforward experimental situations involving only the von Neumann-type experiments (with no 
^ ■ sequential measurements, etc). 

In the early nineties, the blueprints for straightforward Bell tests involving higher dimensional systems were given 
(for a summary see [||). The idea was to use unbiased multiport beam splitters to define the local observables. 
Surprisingly, it turned out, that such observables suffice to reveal the fact that pair of entangled higher dimensional 
^ ( - systems violate local realism more strongly than qubits ||. This result was obtained numerically by employing 
the linear optimization procedures to search for underlying local realistic joint probability distribution that would 
reproduce the quantum prediction (with some noise admixture). The results were confirmed analytically in Refs. ficfl 
and jnj. Later in Ref. |l2|], it was shown that in the case of pairs of entangled higher-dimensional systems, violations 
of local realism are even stronger for non-maximally entangled states. In a parallel research, it has been shown that 
| higher dimensional systems can lead to the GHZ-like paradox without inequalities |lj| , |l4| . 

In view of all these facts, it is tempting to test the strength of violation of local realism by triples of higher 
dimensional systems (starting of course with three qutrits), and that for non-maximally entangled states. 

Since Bell-type inequalities for three qutrit systems are unknown at the moment, it is necessary to invoke the 
numerical algorithm first presented in p5[ . As we shall see, some surprising results can be obtained in this way. 
We show below the result of our numerical analysis. It turns out that 

1. There is a strong violation of local realism (for the standard von-Neumann type measurements) for three qutrit 
systems in the maximally entangled state, however it is not as strong as in the case of the three entangled qubits. 

2. Allowing non-maximally entangled states, the situation changes. We find the three qutrit state which reveals 
correlations much much more resistant to noise, than those for entangled three qubits (maximally entangled 
three qubit states give maximal violation of local realism |l6j, fll7| , rfl8|)- 

In our numerical analysis, we consider a class of pure states of three qutrits in the form of 

10) = £ d gij \g)m) (1) 
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with real coefficients d g ij. The kets \g),\i),\j) denote the orthonormal basis states for the first, second and the 
third qutrit respectively. Three spatially separated observers, Alice, Bob and Cecil, are allowed to perform the 
measurement of two alternative local noncommuting trichotomic observables on the state \ip). We assume that they 
measure observables defined by unbiased symmetric three-port beamsplitters In such a situation the kets in (|l|) 
represent spatial beams, in which the particles can propagate. The observers select the specific local observables by 
setting appropriate phase shifts in the beams leading to the entry ports of the beamsplitters. The overall unitary 
transformation performed by such a device is given by 

1 i2n 

U 3'j = -/g c M—fj) exp(^ ), (2) 

where j denotes an input beam to the device, and j' an output one, and 4>j are the three phases that can be set by the 
local observer (for a more detailed description see ]§| ) . Please note, that the actual physics of the device is irrelevant 
for our theoretical discussion here, thus it suffices just to assume that the observers perform their von Neumann 
measurements in the basis which is related to the "computational" basis of the initial state (|l|) by the transformation 
(||). It is interesting that the unitary transformation for all phase settings leads to a new basis for the local qutrit, 
which is unbiased with the respect to the "computational" one. 

Let us denote Alice's local unitary transformations associated with her device by Ua{<Po), Ua(<Pi), Bob's by 
Ub(xo),Ub(xi) and Cecil's by Uc{S ), Uc($i), where the three component vectors 4>k,Xi,S m (k,l,m = 0,1) de- 
note the set of the phases defining the appropriate observables. The measurement of each observable can yield 
three possible results which we denote by a for Alice, b for Bob and c for Cecil (a,b,c — 1,2,3). The probability 
PQAt{o-k,bi, c m ), that Alice, Bob and Cecil obtain the specific results after performing the unitary transformations 
UA(4>k), Ub{xi) and Uc(S m )i respectively, is given by the following formula 

PQM{a>k,k,c m ) = \{ak\(bi\{c m \UA($k)U B (xi)Uc{S m )\ip)\ 2 

~ 27 ^ 27 Z^g'i'j'^gij *V« 3 ,Ul ff»J 

x cos(f (a k (g - g>) + b t (i - i>) + c m (j - j')) + <fi - 4 + X \- X \ + S m - S m ) , (3) 

where, for instance, <p 9 k denotes the <?-th component of 4>k- 

In the presence of random noise, in order to describe the system one has to introduce the mixed state pf = 
(1 — F)\ip)(ip\ + Fp no i se , where p no i se — ^fl, and I is the identity operator. The non-negative parameter F specifies 
the amount of noise present in the system. In such a case, the quantum probabilities read 



PQM( a k,bi,c m ) = (1 - F)PQ M (a k ,bi,c m ) + 



F 
27 



The hypothesis of local realism assumes that there exists some joint probability distribution Plr{_o-o-, Oi! frcb b\ \ Cq, C\) 
that returns quantum probabilities PQ M (ak,bi,c m ) as marginals, e.g., 

-Pqm (°0j b , c ) 

3 3 3 

= ^2 ^2 ^ -Plh(qq, Qi; bp, bi; c , a). (4) 

ax=l bi—X ci=l 

Please note, that a concise notation of the full set of such conditions can be given by 

PQAl( a k, bl,C m ) 

3 3 3 

= X! ^2 P LR( a o, a i,bo,bi;c ,Ci). (5) 

Qfc + l=l &! + l = l C m + l = l 

where fc+l,Z + l,m+l are understood as modulo 2. For each pure state \ip), one can find the threshold F t h r (the 
minimal value of F) above which such a joint probability distribution satisfying (@) exists (obviously, for any separable 
state F t hr = 0, however this may hold also for some non-separable states). 

There is a well defined mathematical procedure called linear programming that allows us to find the threshold F t h r 
for the given state \ip) and for the given set of observables. We should stress that F t hr found in this way gives us 
sufficient and necessary conditions for violation of local realism. The procedure works as follows. 
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The computation of the threshold F t hr is equivalent to finding the joint probability distribution 
Plr(cio: CLi;bo, bi; cq, c%), i.e., the set of 3 6 of positive numbers summing up to one and fulfilling 8 x 27 — 216 conditions 
given by (||) such that F is minimal. Therefore, F and Plr(clo, Oli bo, b\, cq, c\) can be treated as variables lying in a 
3 6 + 1-dimensional real space. The set of linear conditions (||) and the condition that < F < 1 defines a convex set 
in this space. 

Next, we define a linear function, whose domain is the convex set defined above so that it returns the number F. 
The task of finding F t hr is then equivalent to the search for the minimum of this function. As the domain of the 
function is very complicated, the procedure can only be done numerically (we have used the numerical procedure 
HOPDM 2.30, see g§). 

It is obvious that the F t hr depends on the observables measured by Alice, Bob and Cecil (which in turn depend 
on the set of phases) as well as on the state (indeed, for some unfortunate choices of observables, or the states 
or both, one can have F t hr = 0). Let us clarify, that the task of the linear optimization procedure is each time to 
find the minimal F, for which the relation (^|) can be satisfied by some positive probabilities on its right hand side. 
However, the left hand side of Eq. (|]) depends on the chosen states and observables, and we are interested in the case 
when getting the local realistic model requires a maximal possible admixture of noise, therefore we search for such 
states and observables, for which the minimal F t h r has the largest possible value. There are two possible interesting 
scenarios. We can fix the state and maximize F t hr over the observables. In this way we find the best violation of 
local realism for this given state. Alternatively, we can maximize F t hr over the coefficients defining the state, as well 
as over the observables. This procedure allows us to find the optimal state, and optimal observables measured on this 
state, which can yield the best possible violation of local realism by the class of pure states with real coefficients (Q) . 
Of course, we do not have to limit ourselves to pure states with real coefficients, nor even to pure states but then in 
these cases the number of parameters over which we have to optimize becomes too large for our computers to handle. 

We have applied the procedure described above for the fixed state which we have chosen to be a maximally 
entangled state, i.e., = -^(|111) + 1 222) + |333). Running the program we have found that the threshold amount 
of noise, that has to be admixed to the maximally entangled state, so that the correlations generated by it, for any 
sets of pairs of local settings of the phases, become describable in a local and realistic way, is F t hr = 0.4. The optimal 
observables form the point of view of violations of local realism, i.e., exactly those for which the noise admixture must 
be maximal to get a local realistic model, are defined by the following sets of phases 4>o — (0, 0, |7r), <j>\ — (0, 0, 0); xo = 
(0, 0, 7r), xi = (0, 0, |7r); So = (0, ^ir, 0),Si = (0, n, 0). We can therefore say, that the violation of local realism in this 
case is stronger than for two maximally entangled qutrits, in which case the threshold amount of noise is only 0.304. 
However, it is weaker than the violation by three entangled qubits, for which the threshold amount of noise is 0.5. 

Naturally, one should check whether one can obtain better violations for non maximally entangled states. Therefore 
we have taken the predictions for ([!]), and used a procedure for the maximalization of F t hr over the parameters d g ij 
as well as the observables. 

We have found that the there exists a non-maximally entangled state, and a certain set of local observables, for 
which one requires F t hr = 0.571 noise admixture for the correlations to have a local realistic description. The 
expansion coefficients of the state are given in the table below, whereas the phases defining the optimal observables 
will not be presented here, as they are not easily interpretable. However, for very close local settings given by: 
(fo = (0, §tt, — |tt), 0x = (0, §7T, 0); xo = (0, ±Ztt, -jgn), Xi = (0, 0, 0); S = (0, tt, §§tt), $i = (0, ^tt, -|tt), there is a 
state for which the threshold noise equals 0.570. 



Basis 


|00Q) 


|001) 


|002) 


|010) 


|011) 


|012) 


|020) 


|021) 


|022) 


Coeff 


+0.186 - 


f 0.076 - 


f 0.230 - 


f 0.218 - 


f 0.046 - 


f0.112 - 


+-0.172 - 


f 0.033 - 


f0.247 



Basis 


|100) 1 101) 1 102) |110) | 111) |112) |120) |121) 1 122) 


Coeff 


+0.216 +0.050 +0.110 +0.160 +0.049 +0.236 +0.204 +0.055 +0.235 




Basis 


|200) |201) |202) |210) |211) |212) |220) |221) |222) 


Coeff 


-0.078 +0.406 -0.029 -0.023 +0.385 +0.035 -0.123 +0.393 -0.128 



In summary, we have shown, that for the maximally entangled state three entangled qutrits violate local realism 
stronger than two entangled qutrits (the threshold amount of noise 0.304, see M). The threshold amount of noise to 
get local realistic correlations is 0.4. This violation is not as strong as for three entangled qubits for which one has 
to admix 50% of noise to make the system describable by local realistic theories. However, we can obtain much a 
stronger violation for the non-maximally entangled states. In this case there exists a non-maximally entangled state 
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(see the table) for which Fthr = 0.57, i.e., we have to add 57% of noise before we enter the region in which the state 
admits local and realistic description. 

We must stress, that although for the state given in the table, the threshold amount of noise F t h r = 0.57 gives the 
necessary and sufficient conditions for the existence of local realism, for the measurement of the observables given 
by unbiased symmetric three-port beamsplitters, it does not mean that with a different choice of observables, or by 
allowing complex coefficients in the state ([!]), one cannot increase Fthr- 

Moreover, it is reasonable to expect, that for four or higher number of entangled qutrits the difference between the 
robustness against noise (i.e., the resistance of quantum correlations to classical description) of maximally entangled 
states and non-maximally entangled ones will still increase. Note that, optimal non-maximally entangled state of two 
qutrits (for which the threshold amount of noise is 0.3139) is around 3% more resistant to noise than the maximally 
entangled one (for which the threshold amount of noise is 0.3038). In the case of three entangled qutrits the difference 
between the threshold amount of noise for non-maximally entangled state (0.571) and for maximally entangled state 
(0.4) is about 40%! 

MZ thanks Nicolas Gisin for discussions on this topic. MZ and DK acknowledge the support of KBN, project No. 
5 P03B 088 20. DK, LCK and CHO would also like to acknowledge the support of A*Star Grant No: 012-104-0040. 
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